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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2471 

UNSTEADY  LAMINAR  BOUNDARY-LAYER  FLOW 
By  Franklin  K.  Moore 


SUMMARY 

The  laminar  compressible  boundary  layer  over  an  insulated  flat 
plate  moving  with  a  time-dependent  velocity  has  been  analyzed  in  detail. 
A  group  of  parameters  ariBe  which,  if  large,  provide  that  the  classical 
"starting  from  rest"  solution  applies,  and,  if  small,  that  the  motion 
is  quasi -steady.  These  parameters  relate  to  the  time  required  for 
temporal  changes  to  diffuse  through  the  boundary  layer.  Deviations 
from  the  quasi-steady  velocity  and  temperature  profiles  have  been 
computed . 

Unsteady  laminar  flows  with  pressure  gradient  and,  probably, 
unsteady  turbulent  boundary  layers  are  governed  by  similar  parameters, 
which  may  be  estimated  in  order  to  provide  a  criterion  as  to  whether 
quasi-steadiness  may  be  assumed  for  a  given  problem.  Hie  case  of  a 
fluctuating  velocity  field  passing  over  a  flat  plate  is  discussed  from 
this  point  of  view. 


INTRODUCTION 

The  unsteady  laminar  boundary- layer  flow  over  various  bodies 
starting  from  rest  has  been  analyzed  by  several  investigators,  the 
first  of  whom  was  Blasius  (see  paragraph  65,  reference  l) .  These 
treatments  have  been  concerned  with  such  matters  as  the  onset  of 
separation  on  an  airfoil  and  the  transient  effects  of  impulsive  start ; 
they  employ  an  expansion  of  some  suitable  quantity  in  powers  of  the 
time  elapsed  since  the  start  of  motion.  Thus,  only  the  earliest  phase 
of  motion  is  considered. 

Unsteady  boundary-layer  flow  for  longer  times  elapsed  since 
starting  and  for  higher  speeds  require  consideration.  For  example, 
the  flight  speed  of  the  usual  rocket  missile  varies  continuously  over 
the  entire  trajectory,  and  thus  the  important  boundary- layer  effects 
of  skin  friction  and  heat  transfer  must,  in  principle,  be  regarded  as 
unsteady  for  the  entire  flight.  Recently,  Kaye  (reference  2)  has 
discussed  the  problem  of  unsteady  heat-transfer  effects  associated 
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vith  missile  flight,  and  has  tacitly  assumed  that  at  high  speeds  the 
boundary  layer  responds  with  no  time  lag  to  changes  in  stream 
velocity;  that  is,  that  the  boundary  layer  at  any  instant  is  that  which 
would  be  associated  with  steady  motion  at  the  stream  conditions  pre¬ 
vailing  at  the  same  instant.  Such  a  boundary  layer  is  hereinafter 
called  quasi-steady. 

Another  type  of  problem  of  unsteady  boundary-layer  flow  for  high 
speed  and  long  elapsed  time  concerns  the  response  of  a  boundary  layer 
to  time-dependent  velocity  or  pressure  fluctuations  in  an  otherwise 
steady  outer  potential  flow. 

Tn  the  present  report  of  research  conducted  at  the  MCA  lewis 
laboratory,  the  case  of  compressible  laminar  flow  over  a  semi -infinite 
flat  plate  in  rectilinear  accelerated  flight  through  still  air  is 
treated  in  detail;  the  flight  speed  is  considered  to  vary  with  time  in 
a  continuous  but  otherwise  arbitrary  manner.  The  solution  to  this 
problem,  which  is  intended  as  an  idealization  of  missile  flight,  will 
be  interpreted  to  provide  insight  into  a  more  general  class  of 
unsteady  boundary- layer  flows,  both  laminar  and  turbulent. 

The  notation  used  in  this  report  is  described  in  appendix  A. 


AM1YSIS  OF  FLOW  OVER  FLAT  PLATE  IN  UNSTEADY  FLIGHT 
Equations  of  Motion 

The  motion  to  be  considered  is  that  of  a  semi “infinite  flat  plate 
moving  in  a  straight  path,  normal  to  its  leading  edge  and  in  its  own 
plane.  Aside  from  the  disturbance  due  to  the  motion  of  the  plate,  the 
medium  (air)  is  at  rest.  The  speed  of  flight  may  vary  with  time.  A 
rectangular  coordinate  system  (fig.  1(a))  is  chosen  to  be  at  rest  in 
the  fluid,  such  that  the  plate  moves  With  a  velocity  U(t)  in  the 
negative  x-direction,  and  y  is  measured  normal  to  the  plate.  The 
origin  of  coordinates  is  taken  to  be  the  leading-edge  location  at 
t  =  0. 

The  Prandtl  boundary-layer  assumptions  are  presumed  to  be  valid 
for  this  flow  and,  in  particular,  the  pressure  is  taken  to  be  constant 
throughout  the  fluid.  The  Prandtl  number  and  specific  heats  are 
assumed  to  be  constant.  Thus,  the  appropriate  equations  of 
motion  of  a  compressible  fluid  are: 
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These  equations  may  he  written  in  a  new  coordinate  system 
(fig.  1(h))  fixed  with  reference  to  the  plate,  with  the  origin  at 
the  leading  edge.  Temporarily  denoting  hy  primes  the  new  coordinates 
and  the  corresponding  Telocity  components,  the  appropriate  trans¬ 
formations  are: 


v*  s  v;  u*  s  u  +  U; 

P* 

x1  H  x  +  U  dt;  y 1  =  y;  t’  s  t 

Jo 

Equations  (1)  become: 


P0  =  constant 


Hereinafter,  the  primes  will  not  he  used,  and  x,  y,  t,  u,  and  y 
are  understood  to  he  measured  relative  to  the  coordinate  system  fixed 
in  the  plate  (fig.  1(h)). 

The  temperature-viscosity  relationship  discussed  hy  Chapman  and 
Euhesin  in  reference  3  is  employed;  that  is,  it  is  assumed  that  the 
variation  of  viscosity  with  temperature  may  he  approximated  hy  the 
relation 
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» =  hs  p'»v'h  (3) 

where  C  is  a  constant  obtained  by  matching  equation  (3)  to  the 
Sutherland  formula  at  (for  example)  the  wall.  This  matching  yields, 
according  to  reference  3, 


+  216° 

+  216u  E  ) 


(4) 


The  equations  of  motion  are  simplified  by  the  use  of  a  transformation 
of  coordinates,  similar  to  that  employed  in  reference  4,  and  of  a 
stream  function,  as  follows: 


Y  s 


7 


(p/p*,)  dy;  X  s  x;  T  s  t 


0 


Po»  blf  & 
us  T37  =  3 Y 


fi*  a  r7 

5'"r+5iJo 


p 

"  "  p  \3x.  3y  3tj 


(5) 


J 


Eolations  (3)  and  (5)  transform,  the  equations  (2)  of  motion  to  yield: 


tyT  +^XY  "^X^YY  =  U'(T)  +  ^YYY 

eT  +  Mx  -  =  ^rueYY  + 

The  appropriate  boundary  conditions  on  y  are 

♦  T(X,-)  =  ty(0,Y)  =  U(T) 

♦  y(X,0)  -  0 

♦  (1,0)  =  0 


(6a) 

(6b) 

(7a) 

(7b) 

(7c) 


The  boundary  conditions  on  6  (X, Y)  must  be  chosen  with  regard  to  the 
particular  problem  under  consideration.  (Subscripts  and  primes  denote 
partial  and  ordinary  differentiation,  respectively.) 
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Solutions  for  l|f,  in  the  Case  of  Constant  Acceleration 

If  the  velocity  is  given  hy  the  relation  U(T)  *=  AT,  time  T 
being  measured  from  rest,  the  equations  (6)  may  be  written  in  similarity 
form,  defining 

if  s  a/JTat3/2  g(z,o 

where 

Z  5  Y/zv/W;  {sx/AT2 

These  definitions  enable  equations  (6a)  and  (7)  to  be  written  as 
follows : 


gz  "  \  Zgzz  "  2^ez{;  +  ^zC  ~  8{8zz  e  1  +  8zzz  ^ 

gz(Z,0)  =  gz(®,C)  =  1;  gz(0,0  *  g(0,£)  *  0  (9) 

Equation  (8)  is  obtained  by  considering  U  =  —  C  to  be  a  constant. 

PqB 

Subject  to  the  limitations  imposed  by  this  assumption,  which  is  dis¬ 
cussed  in  appendix  B,  this  analysis  may  be  considered  applicable  to 
compressible  flow. 


Case  £>>1.  -  If  £  s=  X/AT2»1,  the  function  g(Z, { )  may  be 
expanded  in  reciprocal  powers  of  £  as  follows: 


g  -  g0(Z)  +  i  gq  (z)  +  -i-  gz(Z)  +  .  .  . 

C  £2 

fib-0*)  =  ij  e1 ' (®)  =  g2'(»)  =  .  .  .  =  o 


g0'(0)  =  0;  g^©)  =  g2‘(0)  »  .  .  .  «  0 
g0(0)  =  0;  gx(0)  =  g2(0)  =  .  .  .  =  0 

When,  in  equation  (8),  successive  reciprocal  powers  of  £  are  equated, 
this  definition  and  set  of  boundary  conditions  yield  a  sequence  of 
differential  equations,  the  first  of  which  is 

go1  -f  W  =  i  +  go,,‘ 
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and  the 


s0* 

«i  =  e2  =  s3  =  *  *  •  “  0 

TMs  approach  provides  a  solution  equivalent  to  that  for  the  initial 
motion  of  a  flat  plate  with  constant  acceleration  presented  in 
paragraph  66^  reference  1  • 

This  analysis  clearly  pertains  to  the  very  earliest  stages  of 
motion  since  the  function  g(Z,£)  is  expanded,  in  effect,  in  powers 
of  time  elapsed  since  the  start  of  motion.  In  this  connection,  it  may 
be  noted  that  the  inertia  terms  in  the  momentum  equation  drop  out  as  a 
consequence  of  assuming  £»1. 

The  boundary  condition  gg(Z,0)  =  1  has  not  been  applied.  Thus, 

no  consideration  is  given  to  the  growth  of  the  boundary  layer  in  the 
X-direction,  and  the  effect  of  the  leading  edge  is  ignored.  The 
quantity  £  is  equivalent  to  the  ratio  of  the  distance  (X)  aft  of 
the  leading  edge  to  the  distance  f  ±  AT^j  that  the  plate  has 

traveled.  (See  fig.  1(a).)  If  £  »1,  it  is  therefore  proper  to 
consider  that  the  effect  of  the  leading  edge  has  not  yet  been  felt  at 
the  station  X. 

Case  £<<1.  -  If,  however,  the  distance  traveled  is  much  larger 
than  distance  aft  of  the  leading  edge,  that  is,  if  £<<  1,  then 
account  must  be  taken  of  the  growth  with  X  of  the  boundary  layer. 

The  appropriate  expansion  of  g(Z, £)  for  small  values  of  £  is 
therefore 

8  '*/r[80 (^l)  *  Csi(4f) + '  82 ^  + '  ’  ] 

subject  to  the  boundary  conditions 

go’W  =  2;  gq'W  =  g2‘W  -  •  •  •  =  0  , 

gQ'(0)  =  0;  g1*(0)  =  g2'(0)  =  •  •  •  -  0 
gQ(0)  =  0}  g1(0)  =  g2(0)  =  ...»  0 

Substitution  of  this  expansion  of  g(Z,£)  into  equation  (8)  yields  a 
sequence  of  differential  equations  the  first  of  which  Is 


solutions  _z/2 

=  .  1  Z2  +  (rtf^Ze"22/4  +  (Z 2  +  2)  W1/2  I 

2  vJq 


e 


dz 
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eoso*'  +  so’"  =  0 


This  differential  equation  and  the  corresponding  set  of  "boundary  condf 
tions  are  precisely  those  which  apply  in  the  case  of  steady  flow  over 
a  semi-infinite  flat  plate.  (See  reference  5.)  In  the  steady  case , 


the  argument  of 
argument  is 


1  _T_ 


In  the  present  analysis,  the 


1  z  _  1  Y  /  nr  IT  /AT  1  Y  /u[tT 

2A/r  2Aywr/VAT2  2  v  “  2  n/xV  » 

Therefore,  at  any  instant  when  {<<1,  the  boundary- layer  flow  for 
constant  acceleration  is  nearly  that  which  would  obtain  in  steady  flow 
at  a  velocity  equal  to  the  instantaneous  velocity  of  the  accelerated 
plate;  as  C-*o,  the  equivalence  becomes  precise.  The  correction  terms 
arising  when  £  is  small  but  not  zero  will  be  considered  in  a  subse¬ 
quent  section. 


Parameters  for  Unsteady  Flat-Plate  Flow 


Before  the  flow  over  a  flat  plate  with  an  arbitrary  velocity-time 
relation  is  considered,  it  is  desirable  to  determine  the  governing 
parameters.  From  the  important  physical  quantities  X,  V,  U,  U*(T), 

U 1  '  (T ) ,  .  .  .  u(n)(T)  ,  .  .  .  the  following  dimensionless  quantities 
may  be  constructed  (in  addition  to  Reynolds  number  and  Mach  number): 


Quantities  (10)  and  Reynolds  number  are  therefore  the  quantities  expected 
to  govern  the  unsteady  incompressible  boundary  layer  on  a  flat  plate  with 
arbitrary  (but  differentiable)  unsteady  flight  speed.  These  quantities, 
in  the  same  form,  also  govern  the  corresponding  compressible  flow,  as 
will  be  shown  in  the  next  section. 


If  U  =  AT,  as  in  the  previous  section,  then  the  quantities  (10) 

jr 

specialize  to  the  single  parameter  £  -  — - .  This  quantity  £  has 

AT^ 

been  noted  to  he  equivalent  to  the  ratio  of  the  distance  aft  of  the 
leading  edge  to  the  distance  traveled  by  the  plate. 

A  more  direct  physical  meaning  may  be  attached  to  the  quan¬ 
tities  (10),  and  hence  to  the  special  case  £,  as  follows: 
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If  some  physical  quantity  were  suddenly  to  change  at  the  boundary 
of  a  fluid  region  of  thickness  8  in  which  conditions  were  initially 
uniform,  this  change  would  diffuse,  with  time,  throughout  the  region) 
uniformity  would  ultimately  be  established  again.  The  time  for 
this  diffusive  process  to  reach  some  prescribed  degree  of  completion 
will,  of  course,  depend  on  the  distance  5  through  which  it  must  act 
and  on  the  coefficient  of  diffusion  V  of  the  fluid  (kinematic 
viscosity  for  dynamical  changes).  In  particular,  one-dimensional 
diffusion  theory  indicates  that  the  time  required  is  proportional  to 

52/P. 

Temporal  changes  at  the  edge  of  a  boundary  layer  will  be  evened 
out,  or  diffused,  through  the  layer  by  this  mechanism.  Accordingly, 
the  time  S2/v  will  be  taken  as  characteristic  of  the  time  required  to 
diffuse  throughout  the  layer  a  unit  proportional  change  occurring  at 
the  outer  edge.  If  the  boundary  layer  is  substantially  the  quasi-steady 
flow  over  a  flat  plate,  62  is  approximately  proportional  to  PX/U,  and 
the  "diffusion  time"  is  proportional  to  x/u. 

Owing  to  the  unsteadiness  of  the  flow,  changes  in  U  appear  at 
the  edge  of  the  boundary  layer.  The  time  required  for  a  unit  propor¬ 
tional  change  in  stream  velocity  to  take  place  cannot  be  given  by  a 
single  quantity,  but,  rather,  will  be  characterized  by  the  following 
set  of  times: 


Accordingly,  the  ratio  of  "diffusion  time"  x/u  to  the  time  required 
for  a  unit  proportional  change  to  be  imposed  is  characterized  by  the 
following  group  of  quantities: 


These  quantities  are  equivalent  to  the  quantities  (10),  which  may 
thus  be  considered  a  measure  of  the  promptness  with  which  the  boundary 
layer  responds  to  impressed  changes.  Clearly,  then,  if  the  quan¬ 
tities  (10)  are  very  small,  quasi-steady  flow  is  to  be  expected.  If 
the  velocity  is  increasing,  the  boundary  layer  becomes  progressively 
thinner  and  hence  responds  more  and  more  quickly  to  changes  at  its 
outer  edge.  This  effect  clearly  tends  to  establish  quasi-steady  flow. 
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Solution  for  Insulated  Plate  and 
Arbitrary  Velocity-Time  Relation 

Solutions  of  the  compressible-flow  equations  (6)  and  (7)  will  be 
Bought  for  the  case  of  a  flat  insulated  plate  traveling  with  a  speed 
which  may  vary  with  time  in  a  differentiable  though  otherwise  arbitrary 
way. 

If  consideration  is  restricted  to  a  stage  of  motion  where  the 
diffusion  time  is  already  rather  small,  the  stream  function  may  be 
defined  as 


where  0  b 


1  Y 


A fc  and  £n  is  a  function  of  X,T,  the  definition  of 


which  will  be  deduced  in  the  course  of  the  analysis.  Similarly, 


=  i  +Zli*k2(T)  ,  .  .  .) 


(The  "recovery  factor"  is  r (0,4q,C ■•#•)•) 

Substitution  of  definitions  (ll)  and  (12)  into  equations  (6)  and 
(7)  yields 


TTT» 

faao  +  -  "8  ^2  +  2 


V  ^nX  +  *  §  £  nT  +  2  ^  «« 


vao  +  Pr  fra  +  \  Pr(foa)2  -  4Pr  ^2  221  r  +  |  221  Or0  + 


W  ^  W 

I  E  r£”{nT  +  I  ^0  E  r£n^  "  5  ^  E  ft„^) 

n=0  n=0  n=0  *  1 


fa(oc,Cn)  =  2;  fo(0,Cn)  »-f(0,Cn)  *  0 
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Since  the  case  of  an  insulated  plate  is  to  be  treated, 
function  r  thus  satisfies  the  boundary  conditions 


0;  the 


ro(°^n)  =  °‘i  rMn)  =  0  (16) 

Equations  (13)  and  (14)  maybe  considered  correct  only  if  the 
quantity  V  may  be  assumed  constant.  For  constant  wall  temperature 
or  for  incompressible  flow,  equation  (4)  indicates  that  C,  and 
hence  V,  would  be  constant.  However,  under  the  condition  of  no  heat 
transfer,  the  wall  temperature,  and  hence  V,  must  depend  to  some 
degree  on  both  X  and  T.  The  circumstances  under  which  this 
dependence  may  be  neglected  are  discussed  in  appendix  B. 

Equations  (13)  and  (14)  may  be  made  self-consistent  (that  is, 
containing  only  functions  of  0  and  the  various  £  n)  by  defining 


X5Cttt 

U4 


(17) 


These  definitions  correspond  to  the  quantities  (10)  previously  derived 
by  dimensional  reasoning.  If  the  quantities  (17)  are  considered  small, 
the  functions  f  and  r  may  be  eaqianded  as  follows: 


f(0,£n)  s  F(o)  +  Cofo(°)  +  £1*1(0)  .+  •.••+  £o2foo(0)  +  •  •  •  + 

Mlf0l(a)  +  •  •  •  (18) 

r(0,(n)  =  R(0)  +  £oro(0)  +  £ir^(0)  +  •  •  •  +  So2rOo(°)  +  •  •  •  + 

^iroi(a)  +  •  •  •  (19) 

Substitution  of  equations  (17),  (18),  and  (19)  into  equations  (13) 
through  (16)  and  collecting  terms  multiplied  by  the  various  powers  and 
products  of  the  £n  yields,  in  part, 

I"'i  +  1^1 1  -  0  (20a) 

f0««*  +  BTq 1 1  -  2F’f0*  +  3F»'f0  =  -4(2  -  F*)  +  aOF1*  (20b) 

"  +  Ffx  1 1  -  4F'f1*  +  5F*  *f]_  =  4f0»  (20c) 
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F*  (a>)  s=  2j 

F'(0)  =  F(0)  «=  0 

(21a) 

VH  =  fn 

•(0)  =  fn(0)  =  0 

(21b) 

B' '  +  Pr  FR'  +  i  Pr  (F  ’  * ) 2  «=  0 

Cm 

(22a) 

r0'  *  +  Pr  Fr0'  -  2Pr  F*r0  =  Pr 

^-F»«f0«»  +  8B  +  20E*  -  3B*f0j 

(22b) 

rx' '  +  Pr  Frx'  -  4Pr  F'r1  *  Pr 

• 

^-F  *  'f  x  *  *  +  4r0  -  5E  *f 

(22c) 

• 

• 

E*(0)  =  E(<*>)  =  0 

(23a) 

rn’(°)  =  rn(®)  =  0 

(23b) 

Equations  (20a) ,  (21a) ,  (22a) ,  and  (23a)  are  the  equations 
appropriate  to  steady  flow.  Accordingly,  equations  (18)  and  (19)  indi¬ 
cate  that  for  small  values  of  the  £n  the  laminar  compressible  hound  ary 

layer  on  an  insulated  flat  plate  is  nearly  quasi -steady,  with  respect 
to  both  velocity  and  temperature  profiles,  provided  that  fn  and  rn 

.  and  their  derivatives  are  of  unit  order  of  magnitude. 

Magnitudes  of  {n>  -  The  system  of  equations  just  derived  will  be 
valid  only  under  circumstances  in  which  the  Cn  are  small,  and  (to  be 
conservative)  can  be  arranged  in  a  decreasing  sequence.  The  first,  and 

usually  dominant  member  of  this  group,  is  Co  =  £21  in  practical 

U2 

situations,  the  requirement  of  small  Co  will  commonly  be  met.  For 

example,  at  a  distance  of  10  feet  aft  of  the  leading  edge  of  a  body 
having  an  instantaneous  acceleration  of  100  feet  per  second  per  second 
(about  three  times  the  acceleration  due  to  gravity),  the  quantity  Co 

will  be  less  than  10  percent  for  a  velocity  greater  than  100  feet  per 
second  and  less  than  1  percent  for  velocities  greater  than  316  feet  per 
second .  For  the  case  of  constant  acceleration  at  100  feet  per  second 
per  second,  £0  win  te  iess  than  1  percent  after  the  first 

3.16  seconds  of  flight. 

For  constant  acceleration,  Co  is  the  only  nonzero  Cn«  For 
*  flight  at  a  velocity  U  proportional  to  T®, 

C,  =  fc"1)  £2.  to  -  fa~1)(m-2)  >3 

k  1  m  *0  9  *2  p  Co 


nr 
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and  so  forth..  For  this  case,  the  Cn  'taken  in  order  clearly  form  a 
decreasing  sequence  if  C0  itself  is  small  and  m  is  not  too  small. 
Since  the  successive  Cn  involve  increasing  powers  of  velocity  in  the 
denominator,  it  will  he  true  that  the  Cn  will  generally  form  a 
diminishing  sequence  for  moderate  or  high  speeds,  provided  that  U(T; 
is  a  differentiable  function.  The  higher  the  speed,  the  more  sharply 
TJ(T)  may  vary  without  violating  this  requirement. 

Solution  of  momentum  equations.  -  The  solution  of  equation  (20a) 
subject  to  boundary  conditions  (21a)  is  available  in  reference  5.  The 
following  features  of  the  function  F(O)  are  important: 

F*'(0)  =  1.328  (24) 


-  2°]  =  -1-721  <25) 

Equations  (20b)  and  (20c)  subject  to  boundary  conditions  (21b)  have  been 
solved  by  direct  numerical  integration.  The  method  used  is  described  in 
appendix  C.  Values  of  f0(O),  fx(0),  and  their  derivatives  are  presented 

in  table  I.  The.  functions  F‘(o),  f0‘(<>),  and  ^*(0)  are  plotted  in 


figure  2. 

The  wall  shear  stress  T^. 


may  be  obtained  in  the 


] 


following  dimensionless  form: 

2 

ox*^  according  “bo  equation  (24)  and  the  results  presented  in  table  ~Ly 

r~2 =  +  <2-655>  -  ^  ir*  •  •  ](26) 


The  leading  term  of  the  right  side  is  the  quasi-steady  solution.  As 
would  be  expected,  positive  acceleration  provides  a  skin  friction  higher 
than  the  quasi-steady  value. 


Displacement  thickness.  -  The  displacement  thickness  8*  may  be 
defined  as  follows: 
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or,  by  use  of  equations  (5)  and  (ll), 


When  equations  (17)  and  (19)  are  used,  the  following  equation  results: 


When  equation  (25)  and  the  information  in  tables  I  and  II  are  applied 
for  air  having  a  Prandtl  number  of  0.72,  this  expression  becomes 


(27) 


the  leading  term  of  which  is  the  quasi-steady  result.  Positive 
acceleration  results  in  a  boundary  layer  which  is  thinner  than  the 
quasi -steady  boundary  layer. 

Solution  of  energy  equations.  -  The  solution  of  equation  (22a)  is 
available  (for  example,  in  reference  3).  The  remaining  equations  (22b) 
and  (22c),  subject  to  boundary  conditions  (23b),  have  been  solved  for 
Pr  =  0.72  by  direct  numerical  integration.  The  method  used  is  described 
in  appendix  C.  As  a  check  on  numerical  procedure,  equation  (22a)  has 
been  solved  by  the  same  method.  Values  of  B(0),  rQ (O),  ^ (<J),  and 
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their  first  derivatives  are  presented  in  table  II.  The  result 
R(0)  =  0.848  agrees  closely  with  the  value  0.845  given  in 
reference  3.  The  functions  R(o) ,  rQ (o ) j  and  r^(o)  are  plotted  in 

figure  3.  The  recovery  factor  for  Pr  =  0.72  is 

r(0)  =  E(0)  +  Ml  r0(0)  +  2^211  ri(0)  +  .  .  . 

TJ2  U3 


or,  from  table  II, 

r(0)  =  (0.848)^1  -  (2.852)  Ml  +  (3.266)  J  (28) 

the  leading  term  of  which  is  the  quasi-steady  solution  (see  reference  3). 
Positive  acceleration  results  in  a  surface  temperature  lower  than  the 
quasi-steady  value . 


UNSTEADY  FLOW  IN  PRESSURE  GRADIEWT 

The  foregoing  analysis  has  shown  that  the  boundary  layer  on  a  flat 
plate  traveling  with  a  speed  which  varies  with  time  becomes  quasi¬ 
steady  for  vanishing  values  of  the  parameters  XU'/U2,  X2U,,AI’3>  and 
so  forth.  When  these  quantities  are  small  but  not  zero,  proportional 
correction  terms  arise. 

When  the  behavior  of  a  boundary  layer  in  unsteady  flow  with 
pressure  gradient  is  to  be  investigated,  it  would  appear  expedient 
first  to  inquire  whether  the  boundary  layer  may  be  considered  quasi¬ 
steady.  If  this  is  not  done,  needless  effort  may  be  expended  on  an 
essentially  simple  problem.  In  order  to  establish  a  criterion  for 
quasi-steady  flow  in  a  pressure  gradient,  knowledge  is  required  of  the 
significant  parameters  arising  in  the  unsteady  equations  of  motion. 


Integrated  Momentum  Equation 

The  incompressible  laminar  boundary^ layer  equations  for  unsteady 

flow, 

Ut  +  UUX  +  VUy  =  Ut  +  UUX  +  VUyy  .(tt) 

ux;  +  Vy  =  0  (30) 
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(31) 


(32) 


where  8  la  considered  a  measure  of  the  thickness  of  the  boundary 
layer,  as  in  the  Pohlhausen  analysis  (see  paragraph  60,  reference  l). 
The  Telocity  profile  is  written 

g  =  h(x,rj);  tj  s  Z  (33) 

The  differential  equation  (29)  may  be  evaluated  at  y  =  0  to  yield 

ut  +  rox  .  -vfujyV 
whence,  hy  us©  of  definition  (33) ^ 

8%  8^U+ 

-Vi(x>0)  =  — +*W- 


hq^  (x,0)  =  -T  s  -(A  +  {■)  (34) 

where 

S^IXjc 

A  =  111  "p ■  "  ( lfPohlhausen  A",  see  paragraph  60, 

reference  l) 

*  =  WU 

Equation  (34)  indicates  that  the  functional  form  of  h  may  be  written 


may  be  integrated  to  yield 


(08*)  +  ^  (U2&**)  +  Wx  8*  =  V  (|i) 


where  x  and  y  are  measured  along  and  normal  to  the  wall, 
respectively  (fig.  1(c)),  and  where 

08  08 


08 

(i.a)dy;  8«.  H  (i  -  h) 

Jo 


§  =  t(r,n) 


(35) 
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whence, 

5*  =  &H*(r);  &**  =  &H**(T)  (36) 

Substituting  equations  (35)  and  (36)  into  equation  (3l)  yields 


rH*  +  —  (U25H**)  =  h-^O)  (37) 

vuS  1 


Equation  (37)  shows  that,  irrespective  of  the  time  variation  of 
U  and  5,  tbe  boundary  layer  will  be  quasi-steady  if  the  following 
quantities  are  each  sufficiently  small: 


5  /VY  s4Uxt  .  SXt 
VTJ  ’  lit  \vJ>  v2  ’  yZy 


(38) 


As  has  been  mentioned  previously,  the  quantity  52/v  may  be  regarded 
as  a  measure  of  the  time  required  to  diffuse  a  change  through  the 
boundary  layer.  Changes  in  time  of  either  TJ  or  5  may  give  rise 
to  effects  to  be  diffused.  Thus,  the  first  and  last  of  the 
quantities  (38)  may  be  regarded  as  the  ratios  of  diffusion  times  to 
the  times  for  unit  changes  of  stream  velocity  to  occur  at  a  given 
value  of  x  because  of  acceleration  and  rate  of  increase  of  accelera¬ 
tion,  respectively.  The  second  member  of  the  group  is  proportional  to 

i  ^  and  thus  relates  diffusion  time  to  time  for  the  occurrence 
P  o  ot 

of  a  unit  change  in  5. 


Relative  to  a  mass  of  fluid  passing  through  the  boundary  layer, 
temporal  changes  may  occur  by  reason  of  variation  of  U  or  8  with  x. 
Measured  in  this  way,  the  time  for  a  unit  change  in  velocity  to  occur 
is  lAJz.  Forming  the  usual  ratio  with  diffusion  time  yields 

ux,  which  is  the  'Tohlhausen  A".  If  this  parameter  vanishes,  the 

boundary  layer  is  that  appropriate  to  the  case  of  zero  pressure  gradient. 
It  is  interesting  to  note  the  analogy  with  unsteady  motion:  "Quasi- 
flat-plate"  flow  will  occur  even  if  Ux  does  not  vanish,  provided  that 

^  Ux  is  small  enough;  the  requirement  is  that  changes  in  U  diffuse 
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w?11  much  more  rapidly  than  they  occur  in  the  stream.  Accordingly, 
the  third  member  of  the  group  (38)  is  analogous  to  the  fourth  but 
depends  on  "both  Eulerian  and  lagrangian  time  variations  of  TJ0 

It  may  further  be  noted  that  the  quantity  T  (equation  (34))  is 
related  to  the  diffusion  through  the  boundary  layer  of  changes  in  shear 

fj”®®  +Ue  i>kStreflvT0lOOlty  variation:  E™ghly,  the  shear  is  propor- 
ional  to  U/8,  and  because  of  Telocity  variation  alone,  its  total  time 
derivative  is  Ut/8  +  UUx/5.  The  ratio  of  diffusion  time  52/u  to  the 

time  for  unit  proportional  change  of  shear  to  occur  is  thus 


ti(Y±  +  wx V 

V  u\  5  T~J  = 


For  the  case  of  zero  pressure  gradient,  the  parameters  - 

a  /V\  „  t „  ■  ' 

\  P  /  Doth  reduce  to  the  previously  discussed  quantity 


and  the 


parameter 


reduces  to 


Therefore }  it  is  to  "be  expected 


m 


will  usually  be  of  the  same  order  and ,  for 


moderate  or  high  speeds,  will  be  much  larger  than  either 


The  foregoing  discussion  pertains  directly  only  to  incompressible 
flow  although  similar  considerations  might  be  expected  to  apply  in 
compressible  flow.  y . 

For  small,  but  nonzero,  values  of  xU'/lT2  and  x2U* ' /u3 
proportional  correction  terms  of  the  same  order  of  magnitude  arise  in 
formulas  (26)  to  (28)  for  shear  stress,  displacement  thickness,  and 
recovery  factor  in  flow  with  zero  pressure  gradient.  Since  xU'/CJ2 

-  2%t  i  > 


^are  special  cases  of  the  quantities  (38),  it  is 


expected  that  for  cases  involving  pressure  gradient,  the  magnitudes 
of  the  quantities  (38)  indicate  quantitatively  the  order  of  magnitude 
of  departure  from  the  quasi-steady  condition; 

In  a  specific  problem,  the  quantities  (38)  may  be  estimated  as 
follows:  Since  a  criterion  for  quasi-steadiness  is  desired,  the 
boundary-layer  thickness  8  may  be  estimated  on  the  basis  of  quasi¬ 
steady  flow.  If  the  resulting  values  of  quantities  (38)  are  of  such  a 
magnitude  as  to  indicate  nearly  quasi-steady  flow,  the  assumption  made 
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in  estimating  5  is  proper.  In  many  cases,  this  estimate  of  8,  and 
particularly  of  d&/dfc,  which  appears  in  the  second  member  of  the 
croup  (38),  can  be  made  only  by  finding  the  complete  quasi-steady 
solution,  a  gamble  justified  in  view  of  the  great  difficulty  in  finding 
a  complete  unsteady  solution  free  of  such  assumptions. 


* 


Oscillating  Stream  Velocity 

As  an  example  to  which  the  foregoing  discussion  applies,  the 
stream  velocity  over  a  flat  plate  is  represented  as 


U  =  +  u^_  sin  2k  k  ^ 


(39) 


which  denotes  a  flow  wherein  a  sinusoidal  velocity  disturbance  of 
amplitude  ux  and  frequency  k  is  carried  downstream  with  the  mean 

velocity  TJ]_;  the  maximum  magnitudes  of  the  quantities  (38)  are,  from 
equation  (39) , 

ls%t|  8k  U1  -d  -  Ul& 

u!  U1 


VTJ 


max 


8% 


'xt 


max 


'tt 


max 


f  8k\2  ui 

^2rt  Ee6  TJij  ^ 


more  precise  estimate  of 


The  quantity  {—■)  cannot  be  estimated  properly  in  advance. 

Tentatively,  it  may  be  assumed  that  it  is  of  the  same  order  as 
8% tMj,  and  then,  if  it  is  concluded  that  the  magnitudes  of  the 
remaining  quantities  indicate  quasi-steadiness,  the  quasi- steady 
solution  may  be  found  and  then  analyzed  to  provide,  as  a  check,  a 

Thus,  if  ^  «  1  and  (1.5  g)'  «  1,  fl»  in 

question  may  he  regarded  as  quasi-steady.  If,  as  a  r0“g 
the  formula  for  boundary-layer  thictness  for  steady  uniform  laminar 
nJZZ  a  flat  plate  is  used,  and  If  it  is  assumed  that  h  -10  cycles 
per  second,  Ux  =  100  feet  per  second,  and  U]/!^  -  0.1,  then  lor 
approximately  the  first  foot  downstream  of  the  leading  edge, 

„  8k  U1  <  n  n, 

Ee8  u7  Ui  0,01 
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<  0.001 


UNSTEADY  TURBULENT  BOUNDARY  LAYER 

The  physical  interpretation  of  the  governing  parameters  in  IhthImbt. 
flow  would  seem  to  have  sufficient  generality  to  permit  extension  to  the 
case  of  the  turbulent  boundary  layer.  The  most  important  Quantities 
vould  be 


V  .» 


5ii8t 

V  5  % 


where  Mb  to  be  interpreted  as  an  average  diffusion  coefficient. 

The  coefficient  of  diffusion  by  turbulent  mixing  may  be  written  2u‘ 
where  2  is  a  '‘mixing  length”  and  u»  is  related  to  the  intensity  of 
turbulent  fluctuations.  If  the  assumption  is  made  that  l  os  8  and 
u'oc  U,  then  the  parameters  (40)  may  be  written 


In  reference  6  it  was  shown  that  the  parameter  6  may  not  be 

considered  to  govern  the  response  of  a  turbulent  boundary  layer  to  an 
adverse  pressure  gradient.  The  evidence  is  inconclusive  for  favorable 
pressure  gradient.  Thus,  the  assumption  made  as  to  the  form  of  the 
turbulent  diffusion  coefficient  is  of  questionable  value  and  only 
qualitative  significance  should  be  attached  to  the  unsteady  flow 
parameters  6U^./U  and  8-^/u .  In  the  absence  of  pressure  gradient, 

5  may  be  taken  to  be  (very  roughly)  proportional  to  x:  thus 


Therefore,  it  would  seem  that  for  zero  pressure  gradient,  if 

Zllf. 

~^2  is  small  enough,  a  boundary  layer  will  tend  to  be  quasi-steady 
whether  it  is  laminar  or  turbulent. 

Actually,  for  zero  pressure  gradient,  the  quantity  xUt/u2  may  be 

obtained  without  assuming  a  specific  form  for  the  coefficient  of 
turbulent  diffusion.  The  boundary— layer  thickness  5  may  be  thought 
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of  as  the  result  of  lateral  diffusion  of  vorticity  from  a  point  source 
(the  leading  edge  of  the  plate)  moving  with  a  velocity  TJ.  (Compare 
17  -reference  1.)  If  this  process  occurs  without  the  com¬ 
plicating  effects  of  pressure  gradient,  one-dimensional  diffusion  theory 

would  indicate  that  &2  -  where'  V  is  either  a  laminar  (molecular) 

or  turbulent  diffusion  coefficient,  and  x/U  is  the  time  elapsed  since  the 
source  of  vorticity  passed  the  station  x.  Thus,  the  "diffusion  time 
&2/V  appearing  in  the  quantities  (39)  would  be  given  roughly  by  x/U, 

and  the  quantity  reduces  directly  to  xUt/u2. 

Of  course,  the  parameter  xUt/U2  yields  a  rather  unsatisfactory 

criterion  for  quasi- steadiness  in  the  turbulent  case  because  no  indica¬ 
tion  has  been  given  of  the  magnitude  of  departure  from  quasi- steadiness 
associated  with  small  nonzero  values  of  this  parameter. 


CONCLUSIONS 

\  The  irnninar-  compressible  boundary- layer  flow  over  a  flat  plate 
moving  with  a  time-dependent  velocity  has  been  analyzed  in  detai  . 
la  the  course  of  this  analysis,  a  group  of  parameters  arise,the 
magnitudes  of  which  determine  the  nature  of  the  flow  unsteadiness.  If 
the  parameters  are  very  large,  the  classical  start ing-from-re st  s 
tion  applies;  if  very  small,  the ' flow  may  be  regarded  as  nearly  juasi- 
steadv*  that  is,  at  any  i  istant  the  motion  is  nearly  that  which  woul 
be  obtained  in  steady  flow  at  the  conditions  prevailing  at  that  instant. 
The  deviations  of  velocity  and  temperature  profiles  (for  an  insulated 
plate)  from  the  quasi-steady  state  have  been  computed .Relative  to 
quasi-steady  flow,  constant  acceleration  results  in  a  thinner  boundary 
layer  with  greater  skin  friction  and  lower  wall  temperature. 

These  parameters  governing  unsteady  flat -plate  flow  may  be  shown 
to  be  special  cases  of  a  group  of  quantities  which  compare  the  time 
required  to  diffuse  a  unit  proportional  change  in  some  stream  quantity 
to  the  wall  and  the  time  for  a  unit  proportional  change  in  that  quantity 
to  appear  in  the  outer  flow.  For  moderate  or  high  speeds  and  rather 
large  accelerations,  the  first  member  of  this  group  is  commonly  very 
small.  The  generalized  unsteady  flow  parameters  appear  explicitly  in 
the  integrated  momentum  equation  for  boundary-layer  flow  with  an 
unsteady  pressure  gradient. 

In  various  problems  of  unsteady  boundary-layer  motion,  for  example, 
in  the  presence  of  oscillations  in  stream  velocity  or  pressure,  these 
parameters  provide  an  advance  criterion  as  to  whether  the  flow  need  be 
considered  essentially  unsteady  or  whether  the  motion  may  be  considered 
quasi-steady. 
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In  the  absence  of  a  pressure  gradient,  the  unsteady  turbulent 
boundary  layer  is  expected  to  be  governed  by  the  same  unsteady  flow 
parameter  as  governs  the  corresponding  laminar  boundary  layer. 

Lewis  Plight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 

Cleveland,  Ohio,  May  25,  1951. 
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b*,  displacement  and  momentum  thicknesses 

£  dimensionless  coordinate  (x/at2) 

T]  dimensionless  coordinate  (y/6) 

Q 

A 

u 

e 

p 
a 

t 

Subscripts: 

w  evaluation  at  vail  (y  =  0) 

evaluation  in  stream  (y— >») 

Subscript  notation  for  partial  differentiation  is  used  when 
convenient.  Primes  are  used  to  denote  ordinary  differentiation. 
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APPENDIX  B 

ASSUMPTION  OF  CONSTANT  V 

Equations  (20)  and  (22)  are  derived  for  constant  V .  Deriva¬ 
tives  of  V  with  respect  to  time  appear  in  the  following  terms  of 
the  differential  equation  (6)  when  definitions  (ll)  and  (.12)  are 
introduced: 


^vrp  =  "o  U'f _  + 


*T  T_i  dMc*2  r-l  2  fU'  '  "*s\l 

Z=~Tr~  +  —  M*  i 


cr0  + 


at  bnT 
n 


The  quantity  -  CT/C  may  he  obtained  from  equation  (4)  • 

Since  the  flat  plate  is  considered  to  be  insulated  and  since  this 
discussion  is  concerned  only  with  orders  of  magnitude,  may  be 

taken  equal  to  the  stream  stagnation  temperature  ^fl  +  Mo^J* 

When  it  is  noted  that  y  for  air  is  about  7/5, 


Up  C«p 

TT-T 


l  +  ilfc2 


216°  R 


,  1„2  ,  .  1  „  2  ,  216°  R 

1+-M.A1+-M.  +— 


U  CT 

i  Moo2 
o 

,  1  „  2  216°  R 

5  6b 

u*  c 

1  +  iM°°2 

1  +  i  M„2  +  216°  R 

5  e. 

G 

If  M*  =  0,  |jr  =0;  and  if  M**  =  oo, 

rt 

1 4-  I,  U  T  —  n  07  -fVn*  M  «  *1 


JJ  Cj 

TJ*  c 


«  1 


=  i.  if  e. 


is  about  430°  R,  |jjT  q  |  =  0.07  for  Moo**  1,  and  0.25  when  Moo®  2. 

The  errors  in  determination  of  fQ  and  Tq  due  to  neglect  of 

c 

are  of  order  -U-  —2  .  Therefore,  for  reasonable  values  of  Sob  and 
U'  T 

for  subsonic  Mach  numbers,  the  dependence  of  u  on  T  may  properly 
be  neglected.  Of  course,  for  high  speeds  when  1?^  might  be  important. 
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the  correction  to  the  quasi-steady  solution  tends  to  he  of  slight 
importance  since  ^magnitude  of  the  correction  is  governed  bf  such 
quantities  as  XU  /u  ,  which  tends  to  be  small  for  large  velocities. 

Of  il)  «enV  deSStSr  (nfaud 

The  quantity  »x  appears  in  all  cases  in  the  factor  1  +  X  Hi 

and  may  be  neglected  if  X  A  »  X  |  «1.  Equation  (4)  combined  with 

equations  (l2)  and  (l9)  yields,  for  small  t  , 


4  =  + 


2  U'C 


V 


If  rQ(0)  is  taken  to  he  of  unit  order,  it  may  he  shown  that 


neglecting  also  involves  an  error  in 
order  1-^-  — 


and 


of 
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APPENDIX  C 
INTEGRATIONS 

Integration  of  Momentum  Equations 

3,)  « d  its  derivatives  is  known  at  five  successive  points,  « 

fifth-degree  polynomial  may  be  matched  to  the  five  too»  and  a  sixth 
LLoun  value  of  f'  "  (c)  .  This  polynomial  may  be  integrated  t 

yield  f0"(c),  f0'C«),  %<“)  at  the  slxth  p0lf  .ln  *!""!  °* 

the  six  values  of  f0  •”(#).  The  condition  that  l0(c)  and  its 
derivatives  must  satisfy  the  differential  equation  (20b)  at  the  sixth 
point  serves  to  determine  at  that  point  the  value  of  f0  {  l> 

hence  the  values  of  f0(«),  fo'(<0,  *nd  V^'  ^ 

solution  for  five  consecutive  points,  the  solution  may  be  extended  to 
the  next  point. 

The  foregoing  numerical  procedure  may  he  used  for  single-point 
boundary- ralue° problems ,  provided  f(c  is  completely  known  for  the 
f  irst  fiveintervals .  In  this  case,  the  first  five  values  are  found 
S  expanding  f0(o>  and  p(c)  in  a  Taylor's  series  about  c  0. 

The  difficulty  presented  hy  the  fact  that  "boundary  conditions  are 
givens  tro^points  V  -  °,  »)  is 

two-point  problem  into  two  single-point  problems;  as  follows,  me 
solution  fQ(cr)  is  written 

f0(«x)  s  Cfo^1)  (a)  +  f0(2)  (<0  (Cli 

where  C  is  a  constant.  The  quantity  fQW  (a)  is  taken  to  be 
the  solution  of  the  single-point  boundary- value  problem 

f0«  "  '  +  K0(l) '  ’  -  ®  ,fo(1> '  +  3F'  ,f0(l)  “  0  1  (C2) 

f W"(o)-i«  f0W'  (o)-f0W  (o)  =  o  _ 
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and  (a),  the  solution  of 

f0t2)"'  +  Ff0(2)''  -  2F'f0(2)'  +  3F"f0(2)  -  -  4(2  -  F‘)  +  2ffF*' 


^2)‘'  (o) 


(2) 


0  C°)  -  f0(2)  (°)  -  0 


(C3) 


Definitions  (Cl)  to  (C3)  provide  that 

^o(a)  =  Cf0W  (cx)  +  fQ^2^  (cr)  will  satisfy  differential  equation  (20h) 
and  will  yield,  at  a  =  0, 


(1) 


V(o)  -  °j  f0'(0)  -  f0(0)  “  0 

(2) 


If,  after  fQ^A^  and  have  been  determined,  C  is  taken  equal 


to  - 


lim  fo(2)'Ca) 


CT  — >  oo 


(i) 


- ,  the  boundary  condition  fo'(oo)  =  0  will  be 

\o) 


satisfied,  and  the  solution  for  f0(cj)  may  be  written 


fo(a) 


lim  fp  (2)'W 

"-»»  (1)',  , 

fn  (<*) 


f-h)  fn\  4-  -F.  (2) 


(a)  +  fQ'  i  (a) 


The  number  of  intervals  for  which  it  is  necessary  to  compute  fn 

(2)  0 


(1) 


and  f^ 
constant  C 


is  primarily  determined  by  the  accuracy  desired  for  the 


The  solution  of  equation  (20c)  is  also  obtained  by  the  foregoing 
method. 


Integration  of  Energy  Equations 

Equations  (22a)  to  (22c)  are  solved  by  the  same  method  as  is 
used  for  the  momentum  equations,  except  that  a  second-degree  (rather 
than  fifth-degree)  polynomial  is  matched  to  three  successive  values  of 
r!I(a).  This  simplification  is  considered  warranted  by  the  fact  that 
the  energy  equations  are  of  lower  order  than  the  momentum  equations. 


r 
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TABLE  I  -  FUNCTIONS  ASSOCIATED  WITH  VELOCITY  PROFILE 


0 

fo 

f  ? 

I0 

f  * 1 

I0 

fl 

V 

fl" 

0 

0 

0 

3.394120 

0 

0 

-1.878370 

.2 

,057747 

.529447 

1.953432 

-.036758 

-.359945 

*1.653824 

,4 

.194689 

, 802508 

.829535 

-.138740 

*.642866 

-1.147047 

,6 

,365840 

,881837 

.013897 

-.286337 

-.813374 

-.557151 

,8 

f 538543 

,827366 

-.512697 

-.456399 

-,869151 

-.018114 

1.0 

. 691664 

,693876 

*,783065 

*.627640 

-.828805 

.396729 

1.2 

.814136 

,528061 

-.845153 

- . 783550 

-.721068 

.653481 

1.4 

,903266 

.365638 

*,760117 

-.913840 

-.577909 

.753609 

1.6 

,962232 

,229349 

-.594847 

-1.014405 

-.428269 

.724592 

1.8 

,995747 

,128992 

*,409600 

-1.086252 

-.293752 

. 610518 

2.0 

1.014474 

,063862 

-.247975 

*1.133778 

-.186499 

.461844 

2,2 

1,023130 

,026808 

-.130397 

-1.162855 

-.109292 

.313684 

2.4 

1.026420 

,008672 

-.057777 

-1,179289 

-.059187 

.193292 

2 . 6 

1.027286 

,001372 

*.019894 

-1.187871 

-.029603 

.108468 

2.8 

1.019088 

*.000727 

*.003696 

-1.192011 

*.013669 

.055563 

3.0 

1.018913 

000829 

.001437 

-1.193852 

-.005818 

.026086 

3.2 

1.018786 

-.000433 

.002081 

-1.194603 

-.002261 

.011299 

3.4 

1.018711 

*.000377 

,001358 

-1.194883 

-,000770 

.004460 

3.6 

1.018660 

-.000175 

,000668 

*1,195151 

-.000210 

.  .001640 

3*8 

1.018637 

-.000082 

.000276 

*1.195168 

0 

.000600 

4.0 

1,018622 

-.000043 

,000088 

*1,195160 

.000090 

.000250 

4.2 

1,018617 

*.000034 

,000014 

-1.195130 

.000120 

.000160 

4.4 

1.018612 

-,000031 

-,000012 

-1.195100 

,000140 

.000140 

I 
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TABLE  II  -  FUNCTIONS  ASSOCIATED  WITH  TEMPERATURE  PROFILE 


0 

R 

R* 

ro 

ro' 

rl 

V 

0 

0.848000 

0 

-2.418820 

0 

2.769860 

0 

.2 

.835307 

-.126789 

-2.382630 

.374156 

2.678376 

-.846230 

.4 

.797491 

-.250398 

-2.267210 

. 779209 

2.455221 

-1.335023 

.6 

. 735902 

-.362717 

-2.073428 

1.147468 

2.160495 

-1.573702 

.8 

'.  653919 

-.452272 

-1.814377 

1.423885 

1.837448 

-1.634264 

1.0 

.557309 

-•.507224 

r- 1.512653 

1.569377 

1.513939 

.  -1.587060 

1.2 

.453863 

-.519942 

-1.196740 

1.567789 

1.207085 

-1.472302 

1.4 

.352151 

-.490678 

-.894813 

1.431645 

.928167 

-1.309995 

1.6 

.259817 

-.428239 

- . 630473 

1.200147 

. 685411 

-1.113120 

1.8 

.182067 

-.347426 

-.417568 

.926319 

.484160 

-.897970 

2.0 

. 121126 

-.262489 

-.259398 

. 660089 

.326110 

- . 684630 

2.2 

.076520 

-.185622 

-.150731 

.435205 

. 208880 

-.492490 

2.4 

.045900 

-.123264 

-.081592 

.265658 

.126810 

-.334450 

2.6 

.026133 

-.077049 

-.040868 

.149928 

.072500 

-.214940 

2.8 

.014108 

-.045402 

-.018722 

.077734 

.038410 

-.131330 

3.0 

.007206 

-.025243 

-.007696 

.036708 

.019510 

-.061160 

3.2 

.003468 

-.013248 

-.002726 

.015507 

. .010000 

-.034920 

3.4 

.001555 

-.006565 

-.000742 

.005627 

.004880 

-.017880 

3.6 

.000631 

-.003073 

-.000088 

.001555 

.002340 

-.008470 

3.8 

.000209 

-.001358 

.000053 

.000132 

.001200 

-.003640 

4.0 

.000027 

-.000567 

.000032 

-.000242 

.000720 

-.001320 

f  Kd0  .  ] 

i 

L. 10974 

00 

l  rQ<icr  « 

t) 

-  2.99236 

OD 

r  rida 
~0 

=  3.22928 
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y 


(a)  Coordinates  fixed  in  fluid  at  rest. 


(c)  Curvilinear  coordinates;  y  =  0  on  wall. 
Figure  1.  -  Coordinate  systems  used  in  analysis. 
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Figure  2.  -  Functions  associated  with  velocity  profile. 
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0  1  2  3 
Figure  3.  -  Functions  associated  with  temperature  profile. 


